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Abstract 

In this paper we establish a version of the Feynman-Kac formula for 
the stochastic heat equation with a multiplicative fractional Brownian 
sheet. We prove the smoothness of the density of the solution, and the 
Holder regularity in the space and time variables. 



1 Introduction 

Consider the following heat equation on IR'' 

' du , 

-^-Au + c{t,x)u ^^^^^ 
u{0,x) = f{x) , 

where / is a bounded measurable function. If c{t, x) is a continuous function of 
{t,x) e [0,00) X M'', then we have the well-known Feynman-Kac formula (see 
[3]) for the solution to above equation 



u{t, x) = E 



fiSne^p (^j\{t-s,B'i)ds 



where Bf = _Bt + a; is a d-dimensional Brownian motion starting from the point 

X. 

In this paper, we shall extend the above Feynman-Kac formula to the heat 
equation with fractional noise 

dt^^ ^^''dtdxi---dxd (1.2) 
^^(O,^) = f{x) , 
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where W{t, x) is a fractional Brownian sheet with Hurst parameters Hq in time 
and {Hi, . . . ,Hd) in space, respectively. The difference between (|l.ll) and (II. 2p 



is that g^g^l is no longer a function of t and x but a generalized (random) 
function. For this equation, we can still formally write down the Feynman-Kac 
formula 



u{t,x) = E'' 



/(B^exp 6{BU~y)W{dr,dy) 



(1.3) 



where denotes the expectation with respect to the Brownian motion i?f , 
and 5 denotes the Dirac delta function. 

The aim of this paper is to justify the above formula (|1.3p . to show that 
the process u{t, x) is a weak solution to Equation (|1.2p and to establish some 
properties of this process. First, we shall show that Vt.x '■— Jq Jg^d S{Bf_j. — 
y)W{dr,dy) is a well-defined random variable. This will be done in Section 2 
using a suitable approximation of the Dirac delta function, assuming that the 
Hurst parameters satisfy 2Ho + Y^'l^^Hi > d+ 1, Ho > i, and Hi > ^ for 
l<i <d. 

After the definition of the random variable Vt^x, the next problem is to 
show its exponential integrability. With the use of the covariance structure of 
the fractional Brownian sheet W{t,x), we show that u{t,x) has exponential 
moments provided 



i?exp 







<oo, (1.4) 



for any A G K. To show that (|1.4p is true we use a method of Le Gall ([5]) 
together with several other techniques. This is done in Section 3. 

Another main point of this paper is to show that u{t, x) defined by (|1.3p is 
the solution to (|1.2p . Instead of following the classical approach based on Ito's 
formula, which seems complicated in our situation, we use again the approxima- 
tion technique together with Malliavin calculus. This is the contents of Section 
4. 

In the above equation (jl.2p the solution and the noise are multiplied using the 
ordinary product. This gives rise to the Stratonovitch integral when we interpret 
the equation in its integral form. However, there is a number of papers where 
Wick product between the solution and the noise is used, which corresponds to 
the Skorohod integral. We shall obtain the Feynman-Kac formula for Skorohod 
type equation in Section 8. 

Feynman-Kac formula gives an explicit form of the solution. This explicit 
form has several consequences. First, in Section 5, we study the smoothness 
of the density of the probability law of u{t, x) (with respect to the Lebesgue 
measure). In Section 6, we use the Feynman-Kac formula to obtain some Holder 
continuity properties of the solution u(t, x) with respect to t and x. 

The above techniques work for Hi>l/2,i~l,2,...,d. From the condition 
2iJo+X]iLi > "^+1 it follows that Hq must be greater than 1/2 and we cannot 
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allow more than one of the Hi , . . . , Hd to be less than or equal to 1/2. Thus 
if we want to remove the condition Hi > 1/2, i — 1,2, ... ,d, we need d — 1. 
We show in Section 7 that ii Hi = ^ and Hq > | then all previous results 
hold except the smoothness of the density. When d = 1, we can also handle 
the case Hq < 1/2, assuming that the process has a regular spacial covariance. 
This has been done in the companion paper [S] using dilferent techniques. We 
would like to mention that this type of Feynman-Kac formula was mentioned 
as a conjecture in the paper [lO] . 

The appendix contains some technical results used along the paper. 



2 Preliminaries 

Fix a vector of Hurst parameters H — {Hq, Hi, . . . , Hd), where Hi e (57 1)- 
Suppose that W — {W{t,x),t > 0,x E M''} is a zero mean Gaussian random 
field with the covariance function 

d 

E{W{t, x)W{s, y)) = Rh, {s, t) \{ Rh, (x,, yi), 

1=1 

where for any H G (0, 1) we denote by Rh{s, t), the covariance function of the 
fractional Brownian motion with Hurst parameter H, that is, 

RHis,t) = \i\tr+\sr-\t-sr). 

In other words, T4^ is a fractional Brownian sheet with Hurst parameters Hq in 
time variable and Hi in space variables, i = 1, . . . , d. 

Denote by £ the linear span of the indicator functions of rectangles of the 
form {s,t] X {x,y] in ]R_|_ x M'^. Consider in £ the inner product defined by 

d 

{IiO,s] X {0,x] , l{0,t] X io,y])n = Rho (s, t) Rh, {xi ,yi). 

In the above formula, if < we assume by convention that I(Q.xi] = ~I{-xi,o]- 
We denote by H the closure of £ with respect to this inner product. The 
mapping W : /(o.t]x(o,x] ^ W{t, x) extends to a hnear isometry between H and 
the Gaussian space spanned by W. We will denote this isometry by 







VF(^) = / / <l){t,x)Widt,dx), 



if S 7i. Notice that if and are functions in £, then 

X / (t){s,x)ij{t,y)\s-t\'^"''"'^W\x,-y^\'^"^-'^dsdtdxdy, 

J 9?. XR2d 



i=l 
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where an — 11^=0 -^'(^-^^ ^ Furthermore, Ti. contams the class of mea- 
surable functions (p on R+ x M'* such that 

^ d 

/ \^{s,x)^{t, y)\\s- TT _ y^\^H,-^dsdtdxdy < C30. 

We will denote by D the derivative operator in the sense of Malliavin calcu- 
lus. That is, if is a smooth and cylindrical random variable of the form 

fiWi^l),...,W{4>n)), 

(pi G H, f C^(M") (/ and all its partial derivatives have polynomial growth), 
then DF is the 7i-valued random variable defined by 

j=i 1 

The operator D is closable form L^(il) into L^{il; Ti.) and we define the Sobolev 
space D^'^ as the closure of the space of smooth and cylindrical random variables 
under the norm 

\\DF\\,^, = ^E{F^)+E{\\DF\\l). 
We denote by S the adjoint of the derivative operator, given by duality formula 

Ei6iu)F) = Ei{DF,u)^), (2.1) 

for any F G D^^^ and any element u G Ti.) in the domain of 6. The operator 

6 is also called the Skorohod integral because in the case of the Brownian motion 
it coincides with an extension of the Ito integral introduced by Skorohod. We 
refer to Nualart [TT] for a detailed account on the Malliavin calculus with respect 
to a Gaussian process. We recall the following formula, which will be used in 
the paper 

FW{(P)^6iF(b) + {DF,<P)^, (2.2) 

for any (p E H and any random variable F in the Sobolev space D^'^. 

Along the paper C will denote a positive constant which may vary from one 
formula to another one. 



3 Definition and exponential integrability of 

JoJ^dSiBf_^-y)Widr,dy) 

For any e > we denote by Pe{x) the d-dimensional heat kernel 

Peix) = (27r£)"^e"^, X e W^. 
On the other hand, for any 6 > we define the function 



Then, ips{t)pi.{x) provides an approximation of the Dirac delta function 5(t,x) 
as e and S tend to zero. We denote by W^'^ the approximation of the fractional 
Brownian sheet W{t, x) defined by 

W''\t,x)^ f f ips{t-s)Peix~y)W{s,y)dsdy. (3.1) 

Fix X e R'^ and t > 0. Suppose that B — {Bt,t > 0} is a d-dimensional 
standard Brownian motion independent of W. We denote by Bf = Bt + x the 
Brownian motion starting at the point x. We are going to define the random 
variable J^^ d{Bf_^ — y)W{dr, dy) by approximating the Dirac delta function 
5{Bf_, - y) by 

Ali{r,y)= f Mt~s-r)Pe{B:-y)ds. (3.2) 
Jo 

We will show that for any e > and S > the function Al'^ belongs to the 
space Ti. almost surely, and the family of random variables 

V:.f= f I Al'^,{r,y)W{dr,dy). (3.3) 

converges in as e and 5 tend to zero. 

Along the paper we denote by E^{<^{B,W)) (resp. by E^{<^{B,W))) the 
expectation of a functional 3'(i3, W) with respect to B (resp. with respect to 
W). We will use E for the composition E^ E^ , and also in case of a random 
variable depending only on B or W. 

Theorem 3.1 Suppose that 2Ho + X^iLi Hi > d + 1 . Then, for any e > and 
5 > 0, Aj'f defined in iS. 2\) belongs to Ti and the family of random variables 
Vt'^ defined m fS^j converges in to a limit denoted by 

Vt,, = f I 6{BU - y)W{dr, dy) . (3.4) 

Jo JR'' 

Conditional to B , Vt^x is a Gaussian random variable with mean and variance 
Var^(Ft,,)=aH f f \r ~ s\^""-^\{\Bl ^ Bl\^"''^ drds . (3.5) 

Jo Jo 

Proof Fix e, e', 5 and 5' > 0. Let us compute the inner product 

(a^^^^A = an f f p,{B: - y)pAB^, - z) 
xips{t — s — u)(ps'{t — r — v) 

d 

x\u- -Q _ z,\^"'-^dydzdudvdsdr. (3.6) 

i=l 
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By lemmas [9 . 2 1 and 1^75] we have the estimate 

d 

X(ps{t - s - u)tps'{t - r - v)\u - up^""^ Y[ IVi - Zt\'^"'^'^dydzdudv 

i=l 

d 

< C\s^r\^"'^-^l[\Bl-Bl\^"'-\ (3.7) 

1=1 

for some constant C > 0. The expectation of this random variable is integrable 
in [0,t]^ because 

f f \s - r\^""-^T\\B\- Blf"''^ dsdr 
Jo Jo 

= f[E\e"'-' f f \s-r\^"'>+^trH,-d-2^^^^ 
Jo Jo 



k{k + 1) 



<oo, (3.8) 



where 



K^2Ho + J2H,-d-l>0. (3.9) 

i=l 

and ^ is a N{0, 1) random variable. The above computations show that the 
condition 2Ho + X)i;=i > d + 1 is necessary for the random variable Vt.x to 
be well defined. 

As a consequence, taking the mathematical expectation with respect to B 
in Equation (|3.6p . letting e = e' and 5 = 6' and using the estimates p.7p and 
([3:81) yields 



2 

< c. 



This implies that almost surely Al'^ belongs to the space H for all e and S > 0. 
Therefore, the random variables V^'^ = W{Al'^) are well defined and we have 



E-E^(y:::vif) = E-(At::,A^^, 

\ I ri 

For any s ^ r and Bs ^ Br, as e, e', 5 and (5' tend to zero, the left-hand side of 
the inequality (|3.7p converges to \s — r\'^Ha-2 J|^^^ _ B].]^"^ ^ . Therefore, 

by dominated convergence theorem we obtain that E^ E^ {V^'^V^ ) converges 
to St as e, e' , 6 and d' tend to zero, where 

* k{k+1) 
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This implies that V^'^ converges in as e and S tend to zero to a Umit denoted 
by Vt^x- Finally, by a similar argument we show (j3.5[) . ■ 

The next result provides the exponential integrability of the random variable 
Vtx defined in 



Theorem 3.2 Suppose that 2Hq + H.i> d+1. Then, for any A G M, we 

have 



EexpiX / 6{Bf_^ - y)W{dr, dy) ] < oo . 
\ Jo Jr'' 

Proof The proof will be done in several steps. 
Step 1 From (|3.5|) we obtain 



(3.10) 



exp I 



f f\s-rr«-^\{\Bl^Bir^-^dsd^ 
Jo Jo J 



and the scaling property of the Brownian motion yields 



Ee 



XVt. 



where /x = ^an^"^^, where k has been defined in (|3.9p . and 

Y= [' [' \s-rf"°-'f\\Bl-Bif"^-'dsdr. 
Jo Jo ,_i 



(3.11) 



(3.12) 



Then, it suffices to show that the random variable Y has exponential moments 
of all orders. 

Step 2 Our approach to prove that i?exp(AF) < oo for any A G M is 
motivated by the method of Le Gall [9]. For k — 1,...,2"^^ we denote 



2fc- 2 2fc- 1 



2" 



2" 



2fc- 1 2k 



2" 2" 



and define 



I \s-r\^"''-^\\\Bl-Bl\^"^-^dsdr. 



The random variables an,k have the following two properties: 

(i) For every n > 1, the variables Qf„,i, . . . , q;„^2'»-i a-i'e independent. 

(ii) a„,fc = 2-"(«+i)ao, where 

»1 »i d 



ao 



/ / is + rf""-'f\\Bl~Bl\''''-'dsdr, 
Jo Jo 

and _B is a standard Brownian motion independent of B. 
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The condition 2Ho + Y^f^i Hi> d+1 implies that Eao < oo and we deduce 
that 



n=l fe=l 



where the series converges in the sense. 
Step 3 For any integer n > 1, we claim that 



Ea^ < E (^C J^' YllBir^-^ds^ , (3.13) 



for some constant C > 0. In fact, we have 

n d 



Ea^ = E [ l[l[{\s,+t,r'>-'\Bl^-Blf^^-')dsdt. (3.14) 



3 

Using the formula 

we obtain for each i = 1, . . . , d, 



1 r°° 



- Bif"'-^ =T{1- Hi)-^ 

. In \ " 

X / Eexp - ^ \Bl^ - El^ IV, n rr^-dr. (3.15) 
For any n, . . . , t„ > and s\,t\, . . . , Sn, tn € (0, 1), we denote 

q, = (e{bibij^) , Q2=(e{bibi)^) . 

\ J «• / nxn \ J «■ / nxn 

We know that 

^«^pf-El^«. --^USj =det(/ + 2Qi+2Q2)-^. (3.16) 
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Substituting into yields 



< 



X / det(/ + 2Qi + 2g2)-^ TTt-^' dr 

^[0,00)" fj^ 

r(l-i7,)"" / det(/ + 2Qi)-^det(/ + 2Q2)"' H""/^''^^ 

J[0,oo)" .1 



< r(i-i/0"" 



/ dct(/ + 2gi)-5 rfrr^'dr 
/ det(/ + 2Q2)-^ n^/^'^^ 



(3.17) 



where in the above first inequality, we have used the estimates 

(/ + 2Q1 + 2Q2) > ^[(/ + 2gi) + (/ + 2Q2)] > (/ + 2gi)^(/ + 2Q2)^. 

Substituting (PTT)) into (|XTil) . and using the inequality (sj+ij)^^«"^ < sf°"^if""\ 
we obtain 



j=i i=i 



i \2Hi-2 



dsdt 



( 



< 



71 d 



J 0,1 " \ 



\2Hi-2 



ds 



Finally, using Holder's inequality with < p < 2 we get 

/ 



[0,1]" \ = i 



H lb 

< E\{\{\Blf"^-^ds 

•^[0,1]" t=ij=i 



|2Hi-2 



This completes the proof of (|3.13p . 

Step 4 For any A > 0, using (|3.13p and Lemma [^751 in the Appendix we obtain 



Ee^"" < Eexp {^^ H \B\\^"^^^d}j < oo, 



(3.18) 



because p < 1. 

Step 5 Define (f{\) = S(e^(""--^"")). By ([XTS|l . (p{X) < oo for aU A e M. 
Since (p'{0) — 0, for every K > we can find a positive constant Ck such that 
for all Xe[0,K] 

^(A) < I + Ca'A'. 

Define ctn.k = a-n.k — E{an,k)- Fix K > and a G (0, k + 1), where k has been 
introduced in p.9|) . Recall that by property (u) in Step 3, a„,fe = 2-"('^+i)ao. 
For every TV > 2 set &Ar = 2KYY=^^ {I - 2-°(J-i)) and set 6i = 2/\. Then by 
Holder's inequality and properties (i) and (ii) of a„_fc, we have for N >2, 



< 



E'exp 6^ ^ ^ an,, 

y n=l k=l 

n=l fe=l , 



i;exp 2'^(^-i)6jv 51 "w.' 



fc=i 



1-2 



a(JV-l) 



-a(N-l) 



< -Eexp 



^ n=l fc=1 / 



(1_„)(W-1) 



Notice that bN2<^-^^-^^+^^^ < 2K. It follows that, 

(^(&Af2'^(^-i)-('^+i)^)2'"°"""' < (i + Ck6^22««-«-i)^'") 

< exp(C2('^+i-2(-+i))^) 

for a constant C independent of N . By induction we get 



,(l-a)(N-l) 



iSexp 



(TV 2""1 \ / ^ \ 

^A' E E < exp C ^ ^ia+l-2i.+ l))n 

n=l fc=l J \ n=2 j 

< exp (C(l - 2°+l-2(K+l))-lj 
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Letting N tend to infinity and using Fatou's lemma, we obtain 

Eexp{b^(Y ~ EY)/2) < oo, 

where boo — 2-ft'n^i(l ^ 2^"'-') > 0. Since K > is arbitrary, we conclude that 
Eexp(XY) < oo for all A e M. This completes the proof, in view of (13. lip . ■ 

4 Feynman-Kac formula 

We recall that is a fractional Brownian sheet on R+ x with Hurst param- 
eters {Hq, Hi, ... , Hd) where Hi e (|, 1) for i = 0, . . . , d. For any £, J > we 
define ^ 

W''\t,x):= f f ipsit~s)p,{x-y)Wids,dy). 



In order to give a notion of solution for the heat equation with fractional noise 
(jl.2p we need the following definition of the Stratonovitch integral. 



Definition 4.1 Given a random field v — {v{t, x),t > 0,x £ R^} such that 

\v{t, x)\dxdt < oo 

almost surely for all T > 0, the Stratonovitch integral J^^ v{t, x)W{dt, dx) is 
defined as the following limit in probability if it exists 

lim / / v{t,x)W^'^{t,x)dxdt. 

We are going to consider the following notion of solution for Equation (|1.2p . 

Definition 4.2 A random field u = {u{t,x),t > 0,x E K.'^} is a weak solution 
to Equation U.2\) if for any C°° function if with compact support on , we 
have 

u[t,x)(p{x)dx = j f{x)ip{x)dx + — I I u{s,x)A(p{x)dxds 
t 

u{s, x)Lp{x)W{ds, dx), 

almost surely, for all t > 0, where the last term is a Stratonovitch stochastic 
integral in the sense of Definition \4-l\ 

The following is the main result of this section. 

Theorem 4.3 Suppose that 2Hq + Yl'i=i > d + 1 and that f is a bounded 
measurable function. Then the process 

u{t, x) = E^ (^f{B't) exp ^* 5{Bt_, - y)W{dr, dy)^ ^ (4.1) 

is a weak solution to Equation hl.2\) . 
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Proof Consider the approximation of the Equation (|1.2p given by the foUowing 
heat equation with a random potential 



dt 2 

From the classical Feynman-Kac formula we know that 

u''^{t, x) = (f{B't) exp ( f W^'^t - s, B^)ds 



(4.2) 



where is a d-dimensional Brownian motion independent of W starting at x. 
By Fubini's theorem we can write 



W'^\t- s,B^)ds 



\Jo 
t 



Mt - s - r)p,{B'i ~ y)W[dr, dy) ds 



(fsit-s- r)peiB'^ - y)ds W{dr, dy) 



where Vf':? is defined ir 





£,5 
t.X 7 

Therefore, 



f-\t,x)=E^ (/(Snexp(y,^f)) 



Step 1 We will prove that for any a: G M"* and any t > 0, we have 
lim E^\u''^\t,x)~u{t,x)\P = Q, 

e,<5iO 

for all p > 2, where u{t, x) is defined in (|4.ip . Notice that 



E^\u''\t, x)--u{t,x)\P 



E' 



(4.3) 



exp 



(y,^i^)-exp(y,.,) 



< ||/||Li?|exp(y,^f)-exp(l^,,,) 



where Vt^x is defined in (|3.4p . Since exp (v^ ;^^ converges to exp (Vt^x) in prob- 
ability by Theorem 13. 11 to show (|4.3p it suffices to prove that for any A e R 



sup exp (^XVt^^'j 



< oo. 



(4.4) 



The estimate (14. 4p follows from 



(1^ . and (|XTUD : 



E exp XVt 



Eexp ( — 



A' 



H 



<i?exp yC 



< oo. 



\r — s 



\2Ho-2 



drds 



(4.5) 
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step 2 Now we prove that u{t, x) is a weak solution to Equation (|1.2p in the 
sense of Definition 14.21 Suppose ip is a smooth function witli compact support. 
We know that, 



1 



u^' {t,x)(p{x)dx = / f{x)(p{x)dx + — / / M^' (t, x)A(p{x)dxds 

+ f f u^-\t,x)ip{x)W'-^^{s,x)dsdx. (4.6) 
Jo Jr'^ 

Therefore, it suffices to prove that 

lim / / u'^'^ {s,x)(p{x)W'^'^ {s,x)dsdx — / / u{s,x)(p{x)W{ds,dx), 

in probability. From (|4.6I1 and (|4.3p it follows that /g J^^ u^'*(s, x)ip{x)W'^'^ {s, x)dsdx 
converges in to the random variable 

G = / u{t,x)(p{x)dx — I f{x)(p{x)dx — '^ j j u{t, x)Aip{x)dxds 

as e and 5 tend to zero. Hence, if 
t 



B^,s= / iu''\s,x) -u{s,x))ip{x)W^'''{s,x)dsdx 
Jo Jr'^ 

converges in to zero, u(s, a:)(^(x) will be Stratonovitch integrable and 

u{s,x)ip{x)W{ds,dx) — G, 



JW 



which completes the proof. In order to show the convergence to zero of B^ g, we 
will express the product {u^'^{s, x)~u{s, x))W'^'^{s, x) as the sum of a divergence 
integral plus a trace term (see (|2.2p ) 

(u^''^(s,x) - u{s,x))W''-\s,x) 

ft r 

{u^'^{s,x) - u{s,x))(ps{s - r)p^{x - z)5Wr,z 

d 

+ {D{u'''\s,x) - u{s,x)), ipsi-s - ■)pe{x - ■))n ■ 

Then we have 

5. A =/"*/" 6Vi5W, 



e,S = / / (Pr]z0^^r.z 
JR'' 
t 



+ / ip{x){D{u''' {s,x) - u{s,x)),ipsis ~ ■)pe{x - ■))Hdsdx 
Jo Jr'' 

Bis + Bis, (4.7) 
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where 

{u^'^{s, x) — u{s, x))ip{x)ips{s — r)p^{x — z)dsdx, 



4.S,S _ / / /„,6,S/ 



IQ JR'' 

and S{cl)^'^) = /J (pf.'j.SWr^z denotes the divergence or the Skorohod integral 
of (f)'^'^ . For the term g we use the following estimate for the Skorohod 
integral 

E[{Bls)'] < EiW'Wn) + Eim^'^r^^^) . (4.8) 
The first term in (14.81) is estimated as follows 



EiW^'rn) = E[iu^''is,x)-uis,x))iu^-'ir,y)~uir,y))] 
Jo Jw Jo Jr'' 

X(p{x)ip{y){ips{s - ■)pe{x - •), 'Psir - ■)Pe{y ~ ■))ndsdxdrdy. (4.9) 
Using lemmas 19.21 and 19.31 we can write 

{ips{s - ■)Pe{x - ■),(ps{r ~ ■)Pe{y - ■)>■« 



= aH \ I ips{s — <j)tp5{r ~ t)\<t — t]"^^^ ^dadr] 
\Jlo,t]^ J 

X / pe{x - z)p^{y - w)W\z^~ w^\'^"'^'^dzdw\ 

V^'' ti J 

d 

< C\s-r\^"«-^l[\x-y\^"'-\ (4.10) 

for some constant C > 0. As a consequence, the integrand on the right-hand 
side of Equation (|4.9p converges to zero as e and S tend to zero for any s, r, x, 
y due to From we get 

sup sup sup £' (u*'''(s, x))^ < 11/11^ sup sup sup i? exp (2V7;^*) < oo. 

e,i5 xeR''0<s<t e,S xGR'' 0<s<t 

(4.11) 

Hence, from (|4.10p and (14. lip we get that the integrand on the right-hand side 
of Equation (g^ is bounded by C\s - rp^o-^HiLi ~ 2/iP^'~^> ^ot some 
constant C > 0. Therefore, by dominated convergence we get that -E(|10^'''|!|^) 
converges to zero as e and S tend to zero. 
On the other hand, we have 

D{u-'\t,x)) = [f{Bt + x) eMVl^)Al 

where A^'* is defined in p.2p . Therefore, 

E{D{u''\t,x)),D{u'' {t,x)))H 
= E'^E^(f{B}+x)f{Bl+x) 



xexp(F,!i*(Bi) + K!f(i?'))(Af:f(Bi),<f (i?2))„). 
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where and are two independent d-dimensional Brownian motions, and 
here denotes the expectation with respect to {B^,B^). Then from the 
previous results it is easy to show that 

hm EiDiu^'^it, x)),D{u''''^'{t, x)))n 

= E[f{Bl+x)fiB^ + x) 

X exp 1 ^ ^ /* /' |. - rn-^ n - Byr'-'dsdr] 



xan f f \s-r\^"'^-^f\\By-By\'"^-^dsdr 

Jo Jo 



(4.12) 



1=1 

This imphes that u'^'^{t,x) converges in the space D^'^ to u{t,x) as (5 J, and 
e J, 0. On the other hand, we also have 

sup sup sup E \\D{u^'^ {s, x))\\'^ < oo. 

e,5 x£R'^0<s<t 

Then 

E\\D4>^''f^^^ = f I f I E{D{u--\s,x)-u{s,x)),D{u-^\r,y)~u{r,y)))n 
Jo Jr'' Jo Jr 

X(p{x)(p{y){ips{s - ■)pe{x - ■), Lps{r - ■)pe{y - ■))ndsdxdrdy 

converges to zero as e and 6 tend to zero. Hence, by (|4.8p ^ converges to 
zero in as e and 5 tend to zero. 

The second summand in the right-hand side of (|4.7p can be written as 

Bis ^ ^{:x){D{u'''^{s,x)-u{s,x)),Lps{s~ ■)p^{x- ■))udsdx 

Jo Jr"^ 

^ix)E^ {fiB^)exp{v::^){Ali,Ms--)Peix--))n)dsdx 

lo Jr'^ 

^{x)E^ {f{B^) exp (K,,) (S{b:_. - ■), Ms - ■)Pe{x - ■))n) dsdx 

,4 

where 



Bl,-B, 



d 

{Al''Ms--)p,{x~-))n = an I I \r ^ v\'""-'T] \y, ~ z,\'"^-' 



Xips{s~ r)p^{B^ ~ y) 

Xips{s — v)p^{x — z)dydyzdrdrdv, 
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and 



an 



),fs{s~ ■)Pe{x- ■))h 
d 



\2Hi-2 



Lemma 19.21 and Lemma 19.31 imply that 



1/35 (r — v)p^{x — y)dydvdr. 



{d{B:_. - .), Ms - ■)Pe{x - ■))n <C r r'""-' J] jB^p^'-^dr, 

•^0 ^=l 



(4.13) 



(4.14) 



for some constant C > 0. Then, from (|4.13p and (|4.14p and from the fact that 
the random variable J^^^^ p^'^^dr is square integrable because of 

Lemma in^U we can apply the dominated convergence theorem and get that ^ 
and ^ converge both in to 



an 



/(Bf)exp(F,,,) / r 



as e and 5 tend to zero. Therefore B^ ^ converges in to zero as e and i5 tend 
to zero. This completes the proof. ■ 

Corollary 4.4 Suppose 2Ho + J2t=i > d + 1. Then the solution u{t, x) has 
finite moments of all orders. Moreover, for any positive integer p, we have 



p 



E{u{t,xy)=Ei Y[f{Bl+x) 



(4.15) 



X exp 



2 ^ 

j,k=l 



J = l 

t nt 



JQ 



\s — r 



2Ho-2-Q|5J,._5M|2H.-2^^^^ 



where . . . , Bp are independent d-dimensional standard Brownian motions. 

Remark 4.5 A formula similar to for the solution to heat equation 

driven by fractional white noise for Hi = • . . = Hd = ^ is discussed in 



5 Regularity of the density 

In this section we shall use the Feynman-Kac formula established in the previous 
section to show that for any t and x, the probability law of the solution u(t, x) 
of Equation (|1.2[) has a smooth density with respect to the Lebesgue measure. 
To this end we shall show that \\Du{t, x)\\.^ has negative moments of all orders. 
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Theorem 5.1 Suppose that 2Ho + J2i=i Hi > d+1. Fix t > and x eW^. 
Assume that for any positive number p, E\f(Bt + a;)|^^ < oo. Then, the law of 
u{t, x) has a smooth density. 

Proof From Theorem 4.3 we can write 

u{t,x)^E^ [f{B^)exp{Vt,,)] . 

The MaUiavin derivative of the solution is given by 

Dr^yu{t, x) = [f{B't) exp (F*,,) S{B^_, - y)] . 

It is not difficult to show that u{t, x) E 11)°^ . Thus, by the general criterion for the 
smoothness of densities (see [H]), it suffices to show that E (^\\Du(t, x)\\ 
oo for any t > and x £ M.'^. We have 



-2p 

n 



< 



\Duit,x)f^ = E"" [fiBl +x)fiB? + x)e^p{Vt,,{B') + Vt,,{B^)) 
x{6{Blf^-y),S{Bt^,-y))n_ 
= a^S^ [fiBl+x)fiB^+x)e^p{VtAB') + Vt^^{B^)) 

12,4 |2Hi-2, 



Jo Jo 



drds 



where B^ and B^ are independent d-dimensional Brownian motions. By Jensen's 
inequality, we have for any p > 0, 

\\Du{t,x)\\^^^ 

< [aHT'E^ [\f{Bl+x)f{Bl+x)\-Pe^^{-p[VtAB') + VtAB^)\) 

d \~ 



Jo 



|2Hn-2 



^ Bl-r-'^drds 



Hence by Holder's inequality, we obtain 

E\\Du{t,x)\\':;^^ 

< [aH)-' {E\f{Bl + x)f{Bf + x)\-PP^)^ 
X {Eexp{~pp2 [Vt.4B^) + Vt,4B^)]))^ 



X E 



hhh : 



"'0 



-PP3\ P3 



\2Ha- 



''-^XWBl^^^Bl^^^^-^d.ds 



where ~" + ^ + ^ = 1- The first factor Ii is finite by the assumption on / and 
Holder ^s inequality. The second factor is finite by Theorem 13.21 Finally, from 
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Jensen's inequality, we have 



< E 



Jo Jo ,_i 



This completes the proof. 



6 Holder continuity of the solution 

In this section, we study the Holder continuity of the solution to the equation 
(|1.2p . The main result of this section is the following theorem. 



Theorem 6.1 Suppose i/iat 2_ffo+X]f=i ^ d+1 and letu{t,x) be the solution 
of Equation iLS^) . Then u{t, x) has a continuous modification such that for 
any p G (O, -j) (where n has been defined in iff.ff)) ). and any compact rectangle 
I C X R'' there exists a positive random variable Kj such that almost surely, 
for any (s, a;), {t, y) £ I we have 

\u{t, V) - u{s,x)\ < Ki{\t - si" + \y- xp"). 

Proof The proof will be done in several steps. 

Step 1 Recall that Vt,x — Jq J^d S{Bf_j. — y)W{dr, dy) denotes the random 
variable introduced in (13.41) and 



u{t,x)=E'' {f{B^)exp{V{t,,)). 
Set V — Vs.x and V — Vt^y. Then we can write 

E'^luis^x) - u{t,y)\P = i?^|i?^(e^ - e^)^ 

< E^[{E^e^'^^''^^'^'>Y^\E^{V-V)Y^^] 

Applying Minkowski's inequality, the equivalence between the norm and the 
LP norm for a Gaussian random variable, and using the exponential integrability 
property p.lOp we obtain 

E^\uis, x) - u{t, y)\P < C[E^{E^{V - Vf)"] ^ 

<Cp[E^E'^\V-V\Y\ (6.1) 
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In a similar way to (|3.5p we can deduce the following formula for the conditional 
variance oi V ~ V 



\r — V 



Jo 



drdv 



Ja Jo 

^2 J \r - v\^"--^ n \B\_,, - Bl_, + X, - y,\^"^-^drdv 
= aHC{s,t,x,y). 



i=l 



(6.2) 



Step 2 Fix 1 < J < c?. Let us estimate C{s,t, x,y) when s = t, and = 
for all i ^ j- We can write 



C{t,t,x,y) = 2 f f\r-vr'l[E{\e"^~^)E{\e''^-'-\z + e"^-')drdv, 
Jo Jo 

(6.3) 



where z 



and ^ is a standard normal variable. Set /3,- — 2Hj + 1 > 2. 



v^^^""^'^'^ ^ ^ 

By Lemma ESI the factor E{\^\'^"'-^ - \z + ^\'^"'-'^) can be bounded by a 
constant if |r — w| < {xj — yjY, and it can be bounded by C\xj — yj\^'\r — v\~^'/'^ 
if |r — s| > {xj — yjY'. In this way we obtain 

C{t, t,x,y)<C [ \r- vl^-^drdv 

J {0<r,v<t,\r-v\<{xj-yjy} 

+C\xj-yjf^ [ Ir-vl^-^-^^/^drdv 

J {0<r,v<t,\r-v\>(xj-yj)^ 

<C\x,-y,\^\ 

So, from (|6.ip we have 

E'^\u{t,x)-u{t,y)\P < C\x, - y,r. (6.4) 

Step 3 Suppose now that s < t, and x — y. Set 6 ~ X^iLi Bi — d. We have 



C(s, x, x) 



|r — ^drdv 
\r - w|2^°"2(|r ~v\^ -\r-v + t- s\^)drdv 
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10 JO 

< / dT 







The first integral is 0{{t — s)"^"*"^), when t — s is small. For the second integral 
we use the change of variable a ^ r ~ v,v ^ t, and we have 

\r - vf"°-^{\r ~vf -\r-v + t- s\^)drdv 

\^"°-^\\a\^ -\a + t-sf\da 

^2Ho-2^^s _^^^^^_^Y)da 

''\-af"»~^ {{-a - t + sY - {-af)da 
t 

+ / ' {-af""-'^\{-aY -{a + t~sY\da 

Js-t 

=^t[A' + B' + C']. 
For the first term in the above decomposition we can write 

(cr'' - {a + l)^)da 



((T* - {a + iy)da 



<c{t-sr, 

because 2Hq + X^^Li ~^ — 3 < —1. Similarly we can get that 

/oo 
a'H"-^{a^ -{a + l)^)da. 



At last, 



So we have 



Jo 



E"^ \u{s, x) - u{t, y)\P < C{t - s)^P. 



(6.5) 



Step 4 Combining Equation 16.41 and Equation 16.51 with the estimates (|6.ip 
and (|6.2p . the result of this theorem now can be concluded from Theorem 1.4.1 
in Kunita [8j if we choose p large enough. ■ 
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7 Case Hq > j, Hi = I and d = 1 



4 

7.1 Preliminaries 

In this case, all the setup is the same as before except that if (f) and ^ are 
functions in then 

S (14^(0) M^(V')) = (</>, ^)H = «if„ 

X / / (p{s,x)ip{t,x)\s - tf""~^dsdtdx, 
Jo Jq Jr 

where aHo = Ho{2Ho — 1). 

7.2 Definition and exponential integrability of 

f;^f^6{Bl,-y)W{dr,dy) 

Similarly we have the following theorem as well. 

Theorem 7.1 Suppose that Hi = 1/2 and Hq > 3/4. Then, for any e > 
and S > 0, A^'^ defined in 13. S^) belongs to Ti. and the family of random variables 
Vt'^ defined m converges in L to a limit denoted by 



Vt,, = / / 5{B^_^ - y)W{dr, dy) . (7.6) 
Jo Jr 

Conditional to B , Vt^x is a Gaussian random variable with mean and variance 

YaT^{Vt.x)^aHo f f \r - sf"°-H {Br - Bs) drds . (7.7) 
Jo JO 



Proof Fix e, e', 6 and 5' > 0. 



= aH„E^( [ [ p,{B: - y)pAB^r " v) 

\J[0,t]^ JR 

X(ps{t ~ s — u)ips>{t ~ r — v)\u — v\^^°~^ dydudvddsdr^ 

= aHo( f E^p,+,,{B,~ Br) 

X(ps{t — s — u)(ps'{t — r — v)\u — v\^^°~^ dudvdsdr^ 

= aHJ I -^{e + e' + \s-r\)-"' 
\J[o,tY v27r 

xtps{t — s — u)ipsi{t — r — v)\u — v\'^^"~'^dudvdsdr | . 
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By Lemma [ 

(e + e' + |s - r\y^^^ x ips(t -s- u)(ps'{t -r- v)\u - vf"°-^ddudv 

[o,t]^ 

<C|s-r|2^«-5/2. 

Then by the dominated convergence theorem, {yt x^^tx^ ) converges 

to 

V27r 7[o,t]2 

as e, e' , S, and S' tend to zero. This imphes that V^'^ converges in as e 



and S tend to zero to a hmit denoted by Vt,x- As e,S go to 0, £' 



w 



converges to right side of Equation (|7.7p almost surely, and because of the above 
argument, the convergence is also in i^, and this implies Equation (j7.7p . ■ 

7.3 Feynman-Kac formula 

By Proposition 3.3 and Theorem 6.2 in [6 , we have the following theorem. 

Theorem 7.2 Suppose that Hi — 1/2 and Hq > 3/4. Then, for any A G M, we 
have ^ 

Eexp (^A^ Jj{Bf_,-y)W{dr,dy)^ < ex. , 
and for any measurable and bounded function f the process 

u(t, x) = E^ (^f{Bf) exp (^^' ^ 5{BU " y)W{dr, dy)^ ^ (7.8) 

is a weak solution to Equation hl.2\l . 

7.4 Holder continuity 

We also have the following theorem, whose proof is similar to that of Theorem 
6.1. 

Theorem 7.3 Suppose Hi = 1/2, Hq > 3/4 and let u{t,x) be the solution of 
Equation Then u{t,x) has a continuous modification such that for any 

p G (0, Hq — 3/4) and any compact rectangle I C M_|_ x M there exists a positive 
random variable Kj such that almost surely, for any {ti, xi), {t2, X2) (z I we 
have 

\u{t2,X2) - u{ti,Xi)\ < Ki{\t2 -ti\P + \X2 - Xi|2f). 

Proof As in the proof of Theorem 6.1, we have 

E^\u{s, x) ~ u{t, y)\P < Cp \e^E^\V - Vl"^ "'"^ 
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where V = ^ 5{B^-r " z)W{dr, dz) and V = 6(3^^^ ~ z)W{dr, dz). If 
s = t, we can write 

E^E'^IV-V]^ = 2 f f \r - v]^""-^ 
Jo Jo 

xE[S{Br - By) - 5{Br - By+x-y)] drdv 

^ , iTH" K/O/^ — iJi EL- 



2tt Jo Jo 



|r - - e '^i^)drdv. 



For any 2p < 7 < 2Ho - 3/2, we have 1 - e ^"F^ < f ^- . Thus 

V2|r-t;| J 

£;^£;^|y - y|2 < C^|x - ypT. Consequently, we have 

E^\u{t, x) - u{t, y)f < C\x - ypf . (7.9) 
On the other hand, if x — y, 



E^E^\V -V\'^ 



t rt 



\r - v\^"°-^/^drdv 



+ / \r-v\^"°-^{\r-v\-^/'^ -\r-v + t-s\-^)drds , 
Jo Jo 

and by a similar computation as step 3 before, we can get 

E^\u{s, x) - u{t, x)\P < C{t - s)(«o-3/4)p^ (7^^o) 
Combining (|7.9p and (|7.10|) we prove the theorem. ■ 

8 Skorohod type equations and chaos expansion 

In this section we consider the following heat equation on 



M(0,x) = /(x). 



The difference between the above equation and Equation \1.2\ is that here we 
use the Wick product o (see [7] , for example) . This equation is studied in 6 in 
the case Hi ~ ■ ■ ■ = H4 = 1/2. As in that paper, we can define the following 
notion of solution. 

Definition 8.1 An adapted random field u — {u{t,x),t > 0,x £ R''} such 
that E(u'^(t, x)) < 00 for all {t,x) is a solution to Equation iS.l]) if for any 
{t,x) e [0,00) X R'', the process {pt-s{x — y)u{s,y)l[Q t]{s), s > 0,y & R''} is 
Skorohod integrable, and the following equation holds 

u{t,x) ^Ptf{x) + [ [ pt-s{x ~y)u{s,y)SWs,y, (8.2) 



JO JM'' 

where pt{x) denotes the heat kernel and ptf{x) = /^d Pt{x — y)f{]j)dy 
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As in the paper the sohition u{t, x) to (|8.ip admits the following Wiener 
chaos expansion 

oo 

= ^ /„(/„(•, t,a;)), (8.3) 

where /„ denotes the multiple stochastic integral with respect to W and /«(■, t, x) 
is a symmetric element in 7i^", defined explicitly as 

/„(si,?/i, . . .,Sn,yn,t,x) = ^ (8.4) 
In the above equation a denotes a permutation of {1, 2, . . . , n} such that < 

Scr(l) < ■ ■ • < So.(„) < t. 

The following is the main result of this section. 

Theorem 8.2 Suppose that 2Ho + X^iLi Hi > d + 1 and that f is a bounded 
measurable function. Then the process 

u(t,x) — E 



f{Bf ) exp ( f [ 5{BU - v)W{dr, dy) 
\ Jo Jw 

(8.5) 



-Ic^hI I \r-sr'>''l[\Bl-Bl\'''^-'drds 
is a weak solution to Equation ([7 



Proof From Theorem 3.2 we obtain that the expectation E^ in Equation (|8.5p 
is well defined. Then, it suffices to show that the random variable u{t, x) has 
the Wiener chaos expansion (|8.3p . This can be easily proved by expanding the 
exponential and then taken the expectation with respect to B. 

Theorem 3,1 implies that that almost surely S{Bf__. — •) is and element of Ti, 
with a norm given by ()3.4p . As a consequence, almost surely with respect to the 
Brownian motion B, we have the following chaos expansion for the exponential 
factor in Equation 



exp / / SiBf_,-y)Widr,dy) 

-\c.H f r|r-.r-^n|i?;-B:|^^'-^drd.Uf;/„(5„), 

where g„ is the symmetric element in 7i®" given by 

gn{si,yi, . . .,Sn,yn,t,x) = ^ S{Bt_ - yi) ■ ■■S{Bf_ - y„) . (8.6) 
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Thus the right hand side of (|8.5p admits the following chaos expansion 

oo ^ 

U{t,x) = y^-lnihn{-,t,x)), (8.7) 

with 

Kit, X) = [f{Bf)S{Bf_,^ - yi) . . . SiBts„ " Vn)] ■ (8.8) 

This can be regarded as a Feynman-Kac formula for the coefficients of chaos 
expansion of the solution of (|8.1|) . To compute the above expectation we shall 
use the following 

E''[fiBmB^-y)\T.] = / pt^4B^-z)f{z)S{z-y)dz 

= pt-s{B'i-y)f{y). (8.9) 

Assume that < So-(i) < • • • < So-(n) < t for some permutation a of {1, 2, • • • , n\. 
Then conditioning with respect to J^t-s^^^ and using the Markov property of 
the Brownian motion we have 

K{t,x) = E^{E^[5{BU„^^^-y.(n)) 

x---<5(Br „ -y.(i))/(Sf)|.F,_,„,J} 



E^ 



Vain)) ■ ■ ■ KBt-s^.^. - ya(l))Ps^<,rJiB^ ) 



Conditioning with respect to ^t-s„(2) and using (|8.9p . we have 

K{t,x) = E''{E''[SiBt_sl^^^-ya(n)) 

x^(Bf-.„<,,-y.(i))p..,,/(Bf-.„<,,)]|^*-M.,} 



xE'' 



^(2) 



= E' 



- ■ ■ ■ ^(^t-^.(2) - y-(2)) 



XPs.(2)-^.(l)(^r-S„(2) - y<T(l))Ps„(l)/(2/^(l)) 

Continuing this way we shall find out that 

hn{t, X) = pt-s,^,^^ix - ?/^(„)) • • •Pv(2)-V(l)(y<T(2) - y<T(l))Ps,(i)/(2/^(l)) 

which is the same as 
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Remark 8.3 The method of this section can be applied to obtain a Feynman- 
Kac formula for the coefficients of the chaos expansion of the solution to Equa- 
tion W. 



with 



U{t,x) = V —In{hn{-,t,x)) , 

^ — ' ni 

f{B^)6{B^_,^-y,)---S{Bf_,^-y^) 

X exp (^^anJ^J^ \r - sf""-^fl {B^ - Bl\^"'~^ drds^ 

(8.10) 

Remark 8.4 We can also consider Equatoin when d = 1, Hi = 1/2 and 

Hq > 3/4. In this case we see easily that the solution u{t, x) admits the following 
chaos expansion 



n=0 



with 



h„{t,x) 



E' 



/(i?nW-s,-2/i)---W-.„-yn) 

exp(^^aHoJ^ ^ \r-s\^''^'-^S{Br-B,)drds^ 



.11) 



9 Appendix 

Lemma 9.1 Suppose < a < 1, e > 0, a; > 0, and that X is a standard normal 
random variable. Then there is a constant C independent of x and t (it may 
depend on a) such that 

E\x + eX\-°' < Cmin(e-", x~°') . 

Proof It is straightforward to check that K — sup^^Q E\z + < oo. Thus 

E\x + = e-'^E\- + X]"" < Ke'°' . (9.12) 

e 
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On the other hand, 

Elx + eX]-"" = [ \x + ey\-°'e-'^dy 

V27r yj{|:r+£a|>f} 

\x + ey\ "e ^ dy ] . 
{k+£al<f} / 

It is easy to see that the first integral is bounded by Cx^" for some constant 
C. The second integral, denoted by B is bounded as follows. 

1 f {'-^)^ If x2 

B = C- l^r^e — ^dz<C- |z|-"e"i^dz 

^ J\z\<§ ^J\A<% 

e 

Thus we have £;|a; + eX\-°' < C|a;|"". Combining this with (HJH), we obtain 
the lemma. ■ 

Lemma 9.2 Suppose a G (0, 1). There exists a constant C > , such that 

sup / pe{xi + yi)pc'{x2 + y2)\yi - y2\^°'dyidy2 < C\xi - X2\^°'- 
Proof We can write 

/ Peixi+yi)p^'ix2 + y2)\yi~y2\^"dyidy2 = E {\eXi ~ xi - s' X2 + X2\^") ■ 
Thus Lemma 19.21 follows directly from Lemma 19.11 ■ 

Lemma 9.3 Suppose a e (0, 1). There exists a constant C > , such that 
sup / / ipsit- si-ri)ips'{t- S2-r2)\ri-r2\^"dridr2 <C\si- 321"°' 

5,(5' Jo Jo 



Proof Since 

1 _x2 1 
PS{X) > PS{X)I^Q^^^{X) = 77^6 > -^ip^{x), 



the lemma follows from Lemma 



Lemma 9.4 Suppose that 2Hq + Hi > d + 1. Let B^, . . . , B'^ he indepen- 
dent one- dimensional Brownian motions. Then we have 



E(j\'""-^l[\Bl\^''^-^ds\ < 
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Proof Wc can write 



t d 
2Ho-2 



2Hi-2\j^i \2Hi-2 



)drds 



Let X be a standard normal random variable. From Lemma 19.11 taking into 
account that 2 — 2Hi < 1, we have when r < s, 

< CE[\Bl\^"'~\s-r)"^-^) 

< Cr"^'\s ~ r)"'-\ (9.13) 

As a consequence, the conclusion of the lemma follows from the fact that 
Jo Jo 

because 2Hq + Yfi^^ Hi - d - 2 > -1 and Y.Li H,-d> -1. ■ 

Lemma 9.5 Let B^,...,B'^ be independent one- dimensional Brownian mo- 

d f 1 \ 

tions. Ifai e (-1, 0), i = 1, . . . d, and ^ > -2, then E ex^ A /g J] |-B*|"*ds < 

1=1 \ i=l J 

oo for all A > 0. 

Proof The proof is based on the method of moments. We can write 

f^l d \ oo^„ ^ n d 



V ^=l / n=l •^[0,l]"fc=i.=l 

OO 



n=l •'I0<si 



<---<s„<l 



\fc=i 



From Lemma 9.1 since at e (—1,0) we obtain 



E 



= E 



\Bl-Bi_^+Bi_X'\r^^_^ 



< Cisk - Sk^iY 



i/2 



where J^t is the filtration generated by the Brownian motion B\ As a conse- 
quence, taking the conditional expectation of 0^=1 I^Ll"' ^^^^ respect to the 

!Fl^ and J^q, we get 



a-fields.f]_^,.Fl_^, 



vfe=i 
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d 

Let a = cti, then we have 



(„1 d \ oo 



/ 



<---<s„<l] 

(r(a/2 + l)r 



Since a > — 2, the integrals on the right side are equal to 7 , , m 

(n + na/2)T[n + na/2) 

and the series converges for any A > 0. ■ 
Lemma 9.6 For any < a < 1 define 

c„(y) = i?(ier"-iy+ei-"), 

where y > and is a standard norm,al random variable. Then 

C^{y) < C mm {I, {y^ + y^-^)), 
for some constant C > 0. 

Proof Notice first that Caiv) < C where C > is a constant, since linij^^oo E\y+ 
^|~" = 0. On the other hand, we can decompose the function Ca{y) as follows 

C„(y) = f {\x\-"-\y + x\-")e-'^dx 



i=i [ {\x\-"-\y + x\-^)e-'^dx 

^TT {x>0}U{x<-y} 

+ [ {\x\-'' -\y + x\-")e-^dx] 

J{-y<x<0} I 



^ (A + B)., 



2tt 

where A and B denote the first and second integrals, respectively, in the last 
second equation. For integral A we can write 

A = / (a;-«-(a; + y)-«)(e-"'/2-e-(^+^)'/2)da; 
Jo 

poo 2 

< / x-''{x + yf-''[{x + yy -x'']ye-'^dx. 
Jo 

Therefore, 

A< x^-^^lix + y)" -x°']ye-'^dx+ / x'^'Hx + y)" - x°']y'^-'^e-'^ dx. 
Jo Jo 
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For the first integral in the above expression we use the estimate (.x + j/)" — a;" < 
ayx°'~^ and for the second we use {x + y)°' — a;" < y". In this way we obtain 

A < Cy\ 

for some constant C > 0. On the other hand, 

B= x-"(e---e —)dx< / x'^'i^x + y)ydx <Cy^-'' , 

Jo Jo 

for some constant C > 0, which completes the proof of the lemma. ■ 

References 

[1] Atar, R., Viens. F., and Zcitouni, O. Robustness of Zakai's equation 
via Feynman-Kac representations. Stochastic analysis, control, optimiza- 
tion and applications, 339-352, Systems Control Found. Appl., Birkhauser 
Boston, Boston, MA, 1999. 

[2] Carmona, R. A. and Viens, F. G. Almost-sure exponential behavior of a 
stochastic Anderson model with continuous space parameter. Stochastics 
Stochastics Rep. 62 (1998), no. 3-4, 251-273. 

[3] Chow, P.-L. Stochastic partial differential equations. Chapman & 
Hall/CRC, Boca Raton, FL, 2007. 

[4] Freidlin, M. Functional integration and partial differential equations. An- 
nals of Mathematics Studies, 109. Princeton University Press, Princeton, 
NJ, 1985. 

[5] Hu, Y., Lu, F., and Nualart, D. Feynman-Kac formula for the heat equa- 
tion driven by fractional noise with Hurst parameter H < 1/2. Preprint, 
2009. 

[6] Hu, Y. and Nualart, D. Stochastic heat equation driven by fractional noise 
and local time. Probab. Theory Related Fields 143 (2009), no. 1-2, 285-328. 

[7] Hu, Y. and Yan, J. A. Wick calculus for nonlinear Gaussian functional. To 
appear in Acta. Appl. Math. 

[8] Kunita, H. Stochastic flows and stochastic differential equations. Cam- 
bridge University Press, 1990. 

[9] Le Gall, J.-F. Exponential moments for the renormalized self-intersection 
local time of planar Brownian motion. Seminaire de Probabilites, XXVIII, 
172-180, Lecture Notes in Math., 1583, Springer, Berhn, 1994. 

[10] Mocioalca, O. and Viens, F. Skorohod integration and stochastic calculus 
beyond the fractional Brownian scale. J. Funct. Anal. 222 (2005), no. 2, 
385-434. 



30 



[11] Nualart, D. The Malliavin calculus and related topics. Second edition. 

Springer, Berlin, 2006. 

[12] Tindel, S. and Viens, F. Almost sure exponential behaviour for a parabolic 
SPDE on a manifold. Stochastic Process. Appl. 100 (2002), 53-74. 



31 



